Gambler's ruin probability - a general formula 



Guy Katriel* 

Department of Mathematics, ORT Braude College, 
Karmiel, Israel 



Abstract 

We derive an explicit formula for the probability of ruin of a gambler 
playing against an infinitely-rich adversary, when the games have payoff 
given by a general integer-valued probability distribution. 



1 Introduction 

We consider the classical Gambler's Ruin Problem. A gambler, who starts with 
an initial wealth M (M £ N), plays a series of games, where in game t the 
gambler's payoff is a random variable X t (independent for different values of t) 
with 

P(X t = k) = Pk, k £ Z, —v < k < oo. 

v > is thus the maximal possible loss, and we assume p_„ ^ 0. The series 
of games proceeds until the gambler's wealth is less than v, in which case the 
gambler must stop playing, and we say that the gambler is ruined. We are inter- 
ested in the probability of ruin P ru in, which depends on the payoff distribution 
{Pfc}fcL-w an d on the initial wealth M. It is well-known that if the expected 
value of X t is non-positive then P rU m = 1, so that we can assume that we are 
dealing with a favorable game 

oo 

E(X t ) = k Pk> 0. (1) 

k— — v 

Another trivial case is when the gambler's initial wealth is less than v, so that 
ruin occurs immediately, so we henceforth assume M > v. 

We will derive a formula for P ru in, which, surprisingly, seems not to be available 
in the literature, except for some very special cases (see discussion below). 

We define the generating function 

oo 

p( z ) = pk 

fc=- 
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z k . 
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Note that, since Y^k=-uP k = •"■> * ne ser i es J>(z) defines a meromorphic function 
in the unit disk \z\ < 1 of the complex plane, with a unique pole, of order i>, at 
z = 0. 

For integers n > 0,r > 0, the complete symmetric polynomial of order r in 
the variables z\, ■ ■ ■ , z n is dehned as the sum of all products of the variables 
zi, ■ ■ ■ , z n of degree r, that is 

n 

$„, r (zi, • • • ,Z n ) = ^ Y[ Z J- 
ij>0, iiH Hn=i"i=l 

Theorem 1 The equation p(z) = 1 /ias v solutions (counting multiplicities) in 
the unit disk \z\ < 1 of the complex plane, which we denote by r\j (1 < j < v). 
The ruin probability is given by 

v n — 1 

Pruin(M) = ^ ®n,M-n+l(Vl, " " " ) Vn) ]J i 1 ~ Vj)- ( 2 ) 
ra=l j=l 

When the roots r\\ , • • • are distinct, one can use the alternative expression 

p™ m (Af)=i>f n 1 ^- (3) 

We note that empty products, as occurs in @ for n = 1, are to be interpreted 
as 1. The two expressions ©j© are algebraically equivalent, but each has 
its advantages. © allows for the case of multiple roots, and shows that the 
expression on the right-hand side of ([3]) is in fact a polynomial in 771, • • • ,?7„. 
([3]) clearly shows that P r uin(M) is is exponentially decaying as a function of 
M (a well-known fact), with the rate of decay determined by the root r\j with 
maximal absolute value. 

The special case of the gambler's ruin problem in which 

p_i > 0, pa > 0, Pk = 0, |fc| > 2, (4) 
is treated in every textbook of probability theory, and it is shown that 

Prvi n (M)=( — ) M (5) 

V pi 1 

(note that the assumption (P) implies p_i < p\). §5§ is obtained as an especially 
simple case of ©, because when © holds we have v = 1, p(z) = p-\z~ x +po + 
Piz, so that 771 = and ^ gives 

Prmn(M) = $l, M (*?l) = Vl' = ( — )" ■ 

\ pi J 

In general, computing P rU m using © requires numerically solving p(z) = 1 in 
the unit disk \z\ < 1. 

The gambler's ruin problem, going back to Pascal and Huygens (see [T] and [2], 
Sec. 7.5 for history), is usually considered together with a version in which the 
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games also stop if the gambler's fortune exceeds a threshold W (the gambler 
wins, or the casino is ruined). Our problem corresponds to W — oo, and is often 
called gambling against an infinitely rich adversary. The result (O is usually 
obtained as a limit of an expression for finite W when W —> oo. 

Feller ( 4 , Sec. XIV.8) and Ethier (0, Sec. 7.2) discuss the problem with an 
upper threshold W in the case of a payoff distribution {pt}^-^! supported 
both from below and from above, as is natural when W < oo. They show that 
the ruin probabilities can be computed by solving a system of linear equations 
of size W — v + 1 for P ru in(M) (y < M < W). It is also possible to express 
Pruin(M) in the form 



where Xj (1 < j < v + /i) are all the roots of p{z) = 1, with the coefficients cij 
determined by a linear system of equations of size v + /i. Here again there is no 
general explicit expression for the solution, so that efforts have been devoted to 
proving upper and lower bounds for P ru im either in terms of the largest root of 
the equation p(z) — 1 (see 0|4]), or more recently in terms of the moments of 
the payoff distribution [5J [5] . 

One import of Theorem Q] is that the problem with W — oo is simpler than 
the case W < oo, in the sense that an explicit formula for the ruin probability 
can be derived. Comparing ([6]) (for the case W < oo) and (jSJ) (for the case 
W = oo), we see that in both formulas we have linear combinations of the 
M-th powers of roots of the equation p{z) = 1, but in ([5]) the coefficients aj 
must be determined by solving a (potentially large) system of linear equations, 
while in everything is given explicitly. Note also that in © only those roots 
which are inside the unit disk are involved. Our derivation does not use ([6]) or 
any other result for the case W < oo - we deal directly with the case ^ = 00. 
Our approach relies on the construction of appropriate generating functions and 
their analysis. 

A previous result which comes closest to the formula ©, in a very special case, 
was obtained by Skala [8 , who considered the case 



and showed that the formula ([3]) holds. Skala used the system of linear equations 
mentioned above for the case W < 00, and a process of going to the limit 
W -> 00. 

Let us remark that in some contexts 'ruin' would be defined as the event in 
which the gambler's fortune drops to or less, rather than below v. Translating 
the results to this case is trivial: on the right-hand sides of ©,([3]), M will be 
replaced by M + v — 1. 

We conclude this introduction with a numerical example. Let us consider a game 
in which participation costs %v and which a prize which is Poisson distributed 
with mean v + e (e > 0). Then p{z) = z" ly e (l ' +e)(1 ^ z) . We take v = 3 and 
e = 0.01. The zeros of p(z) — 1 in the open unit disk (found numerically using 




(6) 



p- v > 0, > 0, p k = 0, k ^ -v, fi, 
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MAPLE) are then given by: z x = 0.993362, z 2 = -.202699 + .220049i, z 3 = z 2 . 
Using (J3J) we obtain 

Pr«i«(3) = 0.9900, P ram (10) = 0.9456, P ruin (50) = 0.7245, 

Pr«m(100) = 0.5193, P rum (200) = 0.2668, P ruin (500) = 0.0361. 

2 Proof of the ruin probability formula 

Defining St to be the gambler's wealth after the t-th game, we have Sq = M, 
and, for all integer t > 

a J St + Xt, St > v . . 

At+1 - \ ft, < S t < v ■ {7 > 

We define the sequence St so that if ruin occurs at some time to, then St = 
St G {Q, 1, • • • , v — 1} for all t > t , or, in other words, {0, 1, ■ • • , v — 1} are 
absorbing states of the Markov process St- Thus the probability of ruin at or 
before time t is given by X)fc=o P(St — k), and the probability of ultimate ruin 
is 



v-l 

Prmn = 1ml V P(S t = A). (8) 
i— >oo * — » 

For later use, we remark that another way to express (0 is 

{St+1 - k) -{ P(St =k) + £f--^ -P(S t = k-l), 0<k<, ■ (9) 
With each of the random variables St we associate its generating function 

oo 

f t (z) = E(z s t)=Y, P ( S t = k > k - ( 10 ) 

k=Q 

Note that fo(z) = z M . Since the coefficients P(St — k) are non- negative and 
their sum is 1 , the series (|10j) converges uniformly in the closed unit disk | z | < 1 
of the complex plane, and thus defines a continuous function in the closed unit 
disk, which is holomorphic in the open unit disk, with 

\z\<l \Mz)\<l. (11) 

Denoting by 7 the vector space of all power series, f(z) — J^hLo c k zk , we define 
the truncation operator T : V — > T, 



T 



k=0 



k=0 



Lemma 1 The sequence of polynomials (of degree v — 1 ) {r[/ t ]}^ converges 
to a polynomial Q{z) 

Q{z)= hmT[/ t ](z), VzGC, (12) 
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0(1). 



(13) 



Proof: We have 



T[ft](z) 



v-l 

E 

k=Q 



P(S t = k)z\ 



(14) 



and since, for each < k < v—1, P{Sk = t) is increasing as a function t bounded 
by 1, we have convergence of each of the coefficients of T[f t )(z) as t — »■ oo, so 
that the limiting degree-(i^ — 1) polynomial exists, and (|T2l holds. In view of 
©, CE2D and (H3]), we have ■ 

We will find P ru in by finding an explicit expression for Q(z) and taking z = 1. 
We note that the coefficient of z k (0 < k < v — 1) in <3(z) is the probability 
that ruin occurs with final fortune k. 

We now derive a recursive formula relating f t +\ to /(. 
Lemma 2 For all t > 0, 



f t+1 (z) = T[f t ] (z) + p{z) f t (z) - T[f t ] (z) 



(15) 



Proof: We show that, for each integer k > 0, the coefficient of z k in the power 
series ft+i(z), is equal to the coefficient of z k in the power series on the right- 
hand of (fTS"!) . Using the notation \z k ]g(z) to refer to the coefficient of z k in a 
power series g{z) 1 we have 



[A{T[ft]{z)) 



[z k ]f t+1 (z) = P(S t+1 = k), 
0, k > v 



P(S t = k), 0<k<v-l ' 



k — V 



[z k ](p(z)[Mz)-T[f t ](z)]) = Y,Pi- p (St = k-l), 



(16) 
(17) 

(18) 



l=-v 



and using ([9]) we see that the right-hand side of (jT6j) is equal to the sum of the 
right-hand sides of ([T7J) and (H]). ■ 



Introducing the linear operator L : V — > V defined by 

L[f}(z) = T[f}(z)+p(z)[f(z)-T[f}(z) 
Lemma [5] tells us that, for all t > Q, 

ft+i = L[f t }. 

For each w £ [0, 1), we define the function 

oo 

t=o 



(19) 



(20) 
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By (fTTj) , the above series converges uniformly in \z\ < 1 for each w £ [0,1), 
and since each f t is a continuous function in the closed unit disk, holomorphic 
in its interior, F w also has these properties. We note also that the functions 
w — > F w (z) {z fixed) are continuous, uniformly for z in compact subsets of the 
unit disk, which in particular implies that the mapping w — > T[F w ](z) from 
[0, 1) to the space of degree-(^ — 1) polynomials is continuous, a fact that will 
be used below. 

Lemma 3 For any z S C we have 

Q{z}= lim (1 - w)T[F w ](z), (21) 

w— >1 — 

where Q(z) is the polynomial defined by 
Proof: Fix z. Set 

c t = T[f t ](z)-T[f t _ 1 ](z), t>\. 

By Lemma [TJ we have 



5> = -T[f ](z) + lim T[f t )(z) = Q(z) - T[f ](z). 

L J t— too 

Therefore Abel's theorem on power series ([7] Theorem 8.2) implies that 

oo 

lim Y^ctw* = Q{z)-T[fo](z). (22) 

W^l— * 4 

For any w £ [0, 1) we have 

oo oo 

(l-w)J2^T[ft}^) = T[f a ]{z) + Y. wt [ T Vt]{z)-T[f t _ 1 ](z) 

t=0 t=l 

oo 

= T[f Q ]{z) + Y, c t wt - ( 23 ) 
t=i 

Taking the limit w — > 1— in ([23]) . using (|22"1) . we get 

oo 

lim (1-w) y^w t T[f t ](z) = Q(z). 

Since 

oo 

TK](zH5> 4 T[/ t ](z), 
t=o 

we have (gTJ). ■ 

In view of LemmaO we now want to find explicit expressions for the polynomials 
T[F w ](z) (w € [0,1)). We need the following result on the solutions of the 
equation 

p(z) = w~ 1 . (24) 



G 



Lemma 4 For any w € (0,1], {21$ has v roots (counted with multiplicities) 
z = Zj{w) (1 < j < v) in the open unit disk \z\ < 1. 

Proof: We first derive an a-priori bound for all solutions. Restricting the 
function p{z) to the interval (0,1], we have, using the fact that pk are non- 
negative, that p(z) a convex function with lim z ^o+ p{ z ) — +oo, p(l) = 1, 
p'(l) = E(X t ) > 0, so that elementary calculus implies that there is a unique 
z* £ (0, 1) with 

p(z*) = 1, p(z) < 1 for z e (z*,l). (25) 

We claim that any solution \z\ < 1 of p(z) = u> _1 , w G (0, 1], satisfies \z\ < z* . 
Indeed if z satisfies p(z) = w' 1 then 

oo oo 

P (\z\) = J2 pM k > | E p^ fe | = = ™ _1 ^ x < 

k——v k——u 

which, by (2SJ), implies |z| < z*. 
Defining 

oo 

h(z) = z v p{z) = y^p fc _^z fc , 

which is continuous on the closed unit disk and holomorphic in the open unit 
disk, we can write the equation p(z) = w^ 1 in the form 

z v = wh(z). (26) 

When w = 0, z = is a zero of multiplicity v of (f!?6"|) . Since, as shown above, 
all solutions of ([25)1 in the unit disk satisfy the \z\ < z* < 1, we conclude by the 
argument principle that ([26]) has ^ solutions for all w £ [0,1]. H 

Lemma 5 For any to e (0, 1), |z| < 1, 

^ 1/ n— 1 

t[f u) ](z) = - Y]$„, M _„ + i(ziH,-- - ,z„H) TT(2-*i(tfl)), (27) 

n=l j=l 

where Zj(w) (1 <j < v) are the solutions of p{z) = w^ 1 in the unit disk \z\ < 1, 
ordered arbitrarily. 

If we further assume that Zj(w) (1 < j < v) are different from each other then 
we can also write 

Proof: Using the linearity of the operator L defined by (fT9]) . and the relation 
(|20f . we have 

OO OO oo 

t=o t=o t=i 
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which we can write as 

F w (z)=wL[F w ](z) + f (z), 
or, using the definition of L, 



F w {z) = wT[F w ]{z) + wp{z) \F w {z) - T[F w ](z)\ + f (z). 
Solving for F w , we have 

w(l-p(z))T[F w ](z) + f (z) 



F w {z) 



1 — p(z)w 



(29) 



Since F w (z) is holomorphic for \z\ < 1, the numerator of the right-hand side of 
(|29p must vanish whenever the denominator does (otherwise we would have a 
pole), that is 

\z\ < 1, l-p(z)w = w(l - p(z))T[F w ](z) + f (z) = 0. 

In other words, denoting by zj(w) (1 < j < v) the solutions of ((24)) (recall 
Lemma 3]) we have 

w{\ - p{zj{vS)))T\F w \(z j {w)) + / (z;M) =0, 1 < j < v, 

and using wp(zj(w)) — 1, we can rewrite this as 

1 



T[F w ](z 3 (w)) 



1 — w 



fo{zj{w)), l<j<v. 



(30) 



We will now temporarily assume that w is chosen so that Zj(w) (1 < j < v) 
are distinct. This holds for all but a discrete set of values of w € [0, 1), since 
if p(z) — w^ 1 has a multiple root z then p' (z) = 0, so that w" 1 is a critical 
value of p, and the set of critical values of a nonconstant holomorphic function 
is discrete. 

T[F w }(z) is a polynomial of degree v — 1, and since (|3U)) prescribes the values 
of this polynomial at the v points Zj(w) (1 < j < i^), we can use Newton's 
interpolation formula [B] to write this polynomial as 

T[F w }(z) = -^—^ MziWr ■ ■ ,zn(w)]]l(z - Zj (w)), (31) 

71 = 1 j = l 

where /[zi(w), • • • , z n (w)] is the n-th divided difference. We now use the fol- 
lowing result ([B], Theorem 1.2.1): when fo(z) — z and < n < M, 

fa[Zi,--- ,Z n ] =$ n ,M-n+l[zi,--' , Z n ]. (32) 

Substituting into (JSIJ) , we obtain (|27| . 

We have proved ([2T| under the assumption that Zj(w) are distinct, but we may 
now approximate any value of w for which some Zj(w) coincide by values of 
w for which Zj(w) are distinct, and use the continuity of both sides of (|2T|) to 
conclude that it is valid for all w. 
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If instead of using Newton's interpolation formula we use Lagrange's formula, 
we obtain (again under the assumption that Zj(w) are distinct) 



T[F„](z) = T ^jr i f (z j (w))l[ 



z - Zi(w) 



Zj(w)-Zi{w) 

Recalling that fo(z) = z M , we get ((27|) . This formula, of course, does not make 
sense when Zj[w) are not distinct. ■ 

We are now ready to obtain the 

Proof of Theorem\]\ Fixing any \z\ < 1, (|2ip and (|27p imply, 

Q(z) = lim {1 - w)T[F w ](z) (33) 

w— >1 — 



= lim — — — y^$ n . M -n+i(zi(w),--- ,z n (w)) Y\(z-Zj(w)). 

w—>l— 1 — W z — ' -LJ. J 

n=l j=l 

By Lemma 01 we have 



_^ v n—1 
lim y^$ nM ^ n+1 (zi(w),---,z n (w))T\(z-Zj(w)) (34) 

w — y 1 — 1 — w A — ' 

n=l j=l 

^ v n — 1 

V ® n ,M-n+l(Vl, ■ ■ ■ >Vn) Y[(z - Vj), 

1 — w ±L 



n=l 3=1 



where rjj (1 < j < 2/) are the roots of p(z) = 1 in the open unit disk. Combining 
([55)) and we have 

n — 1 

QO) = "j J] $«,M-n+l(?7l, ■•• ,Vn) T[( z - Vj), 

1 — W 

n=l j'=l 

and setting z = 1 and using (fT5|) gives @. 

If we assume that J?i, ■ • • , f? y are distinct, which implies that z\{w), • • • , z„(w) 
are distinct for w sufficiently close to 1, then the same argument, using 
instead of {57}, leads to 



and substituting z = 1 gives (J5J) 
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